Precision measurement of the top quark mass from dilepton events at CDF II 
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We report a measurement of the top quark mass, M t , in the dilepton decay channel of tt — ► 
bl' + vtr~b£~T>t using an integrated luminosity of 1.0 fb _1 of pp collisions collected with the CDF II 
detector. We apply a method that convolutes a leading-order matrix element with detector resolution 
functions to form event-by-event likelihoods; we have enhanced the leading-order description to 
describe the effects of initial-state radiation. The joint likelihood is the product of the likelihoods 
from 78 candidate events in this sample, which yields a measurement of M t = 164.5 ± 3.9(stat.) ± 
3.9(syst.) GeV/c 2 , the most precise measurement of Mt in the dilepton channel. 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk, 12.15.Ff 



The top quark, the weak isospin partner of the bottom 
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quark, is the most massive of the known fundamental 
particles. The top quark mass, M t , is a fundamental pa- 
rameter in the standard model. Precise measurements of 
M t along with those of other standard model parameters 
can be used to place constraints on the mass of the Higgs 
boson ll| and on particles in extensions to the standard 
model [2|. Currently, the Tevatron collider at Fermilab 
is the only accelerator capable of producing top quarks, 
where they are primarily produced in pairs. The dilepton 
channel, including decays with two charged leptons in the 
final-state (ti — > W + bW~b — > bl' + viiM~Vi), has a small 
branching fraction but has the fewest jets in the final- 
state, giving a smaller dependence on the calibration of 
the jet energy scale and less ambiguity in jet-quark as- 
signments. Nevertheless, discrepancies between measure- 
ments in different decay channels could indicate contribu- 
tions from physics beyond the standard model [3j . Previ- 
ous measurements of M t in the dilepton channel 0, @, @] , 
while statistically limited, have yielded lower values than 
measurements in other decay channels 0, S, 0, E3] • 

The dilepton channel poses unique challenges in re- 
constructing the kinematics of ti events as two neutrinos 
from W decays escape undetected. Measurements of M t 
in this channel made using Run I data [H, H and re- 
cent measurements made using Run II data utilize 
methods that make a series of kinematic assumptions 
and integrate over the remaining unconstrained quan- 
tities. The greatest statistical precision, however, was 
achieved throu gh the application of a matrix-element 
method [II, HH which makes minimal kinematic as- 
sumptions, instead integrating the leading-order matrix- 
element for ti production and decay over all uncon- 
strained quantities. The first application of this method 
to the dilepton channel by the CDF collaboration [i , 14 1 
used 340 pb _1 of Run II data. 

This Letter reports a measurement using an enhanced 
version 11511 of the matrix-element method described in 



Ref. [lj] . The enhanced method accounts for initial-state 
radiation from the incoming partons and has substan- 
tially improved statistical power. This measurement uses 
data collected by the CDF II detector between March 
2002 and March 2006 corresponding to an integrated lu- 
minosity of 1.0 fb^ 1 and includes the 340 pb _1 used in 
Ref. Q. 



The CDF II detector 16] is a general-purpose detec- 
tor, designed to study pp collisions at the Tevatron col- 
lider. The charged particle tracking system consists of a 
silicon microstrip tracker and a drift chamber, both im- 



mersed in a 1.4 T magnetic field. Electromagnetic and 
hadronic calorimeters surround the tracking system and 
measure particle energies. Drift chambers located out- 
side the calorimeters detect muons. 

The data used in this measurement are collected using 
the same triggers as in Ref. [14| . After events passing 
the trigger requirement are reconstructed, we impose the 
selection criteria defined as "DIL" in Ref. [l7| to isolate 
the dilepton candidates. These selection cuts yield 78 
candidate events. 

We express the probability density for ti decays as 
P s (x|M t ), where M t is the top quark mass and x rep- 
resents the lepton energy, jet energy, and $t [IH mea- 
surements. We calculate P s (x\M t ) using the theoretical 
description of the tt production and decay process with 
respect to x; P s (x.\M t ) is proportional to the differential 
cross section, da(M t )/dx. 

We evaluate P s (x\M t ) by integrating over quantities 
that are not directly measured by the detector, such as 
neutrino momenta and quark energies. While quark ener- 
gies cannot be directly measured, they can be estimated 
from measured jet energies. We integrate over quark en- 
ergies using a parameterized transfer function Wj e t(p, j), 
which is the probability of measuring jet energy j, given 
quark energy p. 

As in Ref. we assume that lepton energies and 

quark angles are perfectly measured, that incoming par- 
tons are massless and have no transverse momentum, and 
that the two highest energy jets in the event correspond 
to the b quarks from ti decay. Unlike in Ref. we do 
not assume zero transverse momentum of the ti system, 
p", which would require no initial-state radiation. In- 
stead, we infer likely values of p" from unclustered trans- 
verse energy [l9| and jets that are not the two most ener- 
getic in the event. We parameterize the relation between 
these measured quantities and p" as a transfer function, 
W PT (pj, , U) , where U is a sum of the unclustered trans- 
verse energy and sub-leading jet transverse energies in 
an event. Both Wj et (p,j) and W PT (p^,U) are estimated 
using ti events generated with HERWIG [2(j and the CDF 
II detector simulation [2l| . This description of the initial 
state radiation improves the expected statistical uncer- 
taint y b y 10% compared to the technique described in 
Ref. [13. 

The effect of the above assumptions on the final mea- 
surement is estimated using Monte Carlo simulation. 
The expression for the probability density at a given mass 
for a specific event can be written as 



P a (x\M t ) = j= f d$\M t t(q hPi ;M t )\ 2 ]J W je t(pk,jk)W PT (pgU)f PDF ( qi )f PDF (q2), (1) 



k=l,2 

r 



where the integral d& is over the eight remaining un- 



constrained momenta of the initial and final-state par- 
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tides, qi and qi are the incoming parton momenta, pi 
are the outgoing lepton and quark momenta, fpnF(<li) 
are the parton distribution functions (PDFs) |22j and 
M. t j(qi,Pi \ M t ) is the leading-order it production and de- 
cay matrix element as defined in Refs. [23L l2~i ] for the 

The term 1/N is 



25] 



process qq — > tt — > bvvgbi . Vy 
defined such that the probability density satisfies the nor- 
malization condition, j dx P s (x\Mt) = 1. The probabil- 
ity for both possible jet-parton assignments is evaluated 
and summed. 

In addition to it production, we calculate the probabil- 
ity for dominant background processes. The final event- 
by-event probability is then P(x.\M t ) = P s (x\M t )p s + 
Pl >1 (x.)pf )1 + Pfc 2 (x)pb 2 • • •, where p s and p\ )i are deter- 
mined from the expected fractions of signal and back- 
ground events (see Table fl} . To determine the , we 
numerically evaluate background matrix elements using 
algorithms adopted from the ALPGEN [26[ generator. We 
calculate probabilities for the following background pro- 
cesses: Z/-y* — > ee,fj,fx plus associated jets, W+ > 3 jets 
where one jet is incorrectly identified as a lepton, and 
WW plus associated jets. We do not calculate prob- 
abilities for Z — > tt or WZ, comprising 11% of the 
expected background. Studies indicate that use of the 
background probabilities improves the expected statisti- 
cal uncertainty by 10%. The posterior probability for the 
sample is the product of the event-by-event probabilities. 
The mean of the posterior probability, P(M t ), is the raw 
measured mass, M t rat0 , and its standard deviation is the 
raw measured statistical uncertainty, AM[™. Both are 
subject to corrections, described below. 



TABLE I: Expected numbers of signal and background events 
for a data sample of integrated luminosity of 1.0 fb _1 . The 
number of expected tt is given for <x t j = 9.1 pb, which corre- 
sponds to Mt = 165 GeV/c 2 . Other backgrounds are negli- 
gible; expected signal and background numbers have an ad- 
ditional correlated uncertainty of 6% from uncertainty in the 
sample luminosity. 



Source 




Events 


Expected it (M t = 165 GeV/c z ) 




63.4 ±1.7 


Expected Background 




26.9 ±4.8 


Drell-Yan (Z/7* -> ££) 


13.1 ± 


4.4 


Misidentified Lepton 


8.7 ± 


1.5 


Diboson (WW/WZ) 


5.1 ± 


1.0 


Total Expected (Mt = 165 GeV/c^) 




90.3 ±5.1 


Run II Observed 




78 



To test the performance of our method we perform 
Monte Carlo experiments of signal and background 
events. Signal events are generated using herwig for top 
quark masses ranging from 155 GeV/c 2 to 195 GeV/c 2 . 
Background events are modeled using observed events 
in the case of background due to misidentified leptons, 
ALPGEN-simulated events in the case of Z/-f* — > ee,/x/z, 
and PYTHiA-simulated [27| events in the case of Z/j* — > 



tt, WW, WZ, and ZZ. The numbers of signal and 
background events in each Monte Carlo experiment are 
Poisson-fluctuated values around the mean values given 
in Table |TJ The estimate for the it signal at varying 
masses is evolved to account for the variation of cross- 
section and acceptance. The response of the method 
for these Monte Carlo experiments is shown in Fig. Q] 
(left). While the response is consistent with a linear de- 
pendence on the top quark mass, its slope is less than 
unity due to the presence in the sample of background 
events for which probabilities are not calculated. Correc- 
tions, M t = 178 GeV/c 2 ± (M[ aw - 176 A GeV/c 2 )/0.83 
and AM t = AM[™ /0.83, are derived from this response 
and applied to values measured in the data. 

The width of the pull distributions in these Monte 
Carlo experiments, shown in Fig. [T] (right), where pull 
is defined as (M t - M* rue )/AM t , indicates that the sta- 
tistical uncertainty is underestimated by a factor of 1.17, 
after applying the corrections described above. This re- 
sults from the simplifying assumptions described above, 
made to ensure the computational tractability of the inte- 
grals in Eq. [TJ The largest effects [lij are the leading two 
jets in an event not resulting from 6-quark hadronization, 
imperfect lepton momentum resolution, imperfect jet an- 
gle resolution, and unmodeled backgrounds. Correcting 
by this factor of 1.17, we estimate the mean statistical 
uncertainty to be 5.0 GeV/c 2 if M t = 175 GeV/c 2 or 
4.2 GeV/c 2 if M t = 165 GeV/c 2 . 
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FIG. 1: Left: Mean measured Mt in Monte Carlo experiments 
of signal and background events at varying top quark mass. 
The solid line is a linear fit to the points. Right: Pull widths 
of Monte Carlo experiments of signal and background events 
at varying top quark mass. The solid line is the average of all 
points, 1.17 ±0.02. 

Applying the method and corrections described above 
to the 78 candidate events observed in the data, we mea- 
sure M t = 164.5 ±3.9(stat.) GeV/c 2 . Figured] shows the 
joint probability density, without systematic uncertainty, 
for the events in our data set. 

The measured statistical uncertainty is consistent with 
the distribution of statistical uncertainties in Monte 
Carlo experiments where signal events with M t = 
165 GeV/c 2 are chosen according to a Poisson distri- 



bution with mean N f - 



63.4 events. This number of 
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FIG. 2: Joint posterior probability density as a function of the 
top quark mass for the 78 observed candidate events, after all 
corrections. Systematic uncertainties are not shown. 



TABLE II: Summary of systematic 


uncertainties. 


Source 


AM t (GeV/c 2 ) 


Jet energy scale 


3.5 


Limited background statistics 


0.7 


PDFs 


0.8 


Generator 


0.9 


Background modeling 


0.2 


FSR modeling 


0.3 


ISR modeling 


0.3 


Response correction 


0.6 


Sample composition uncertainty 


0.7 


Lepton energy scale 


0.1 


Total 


3.9 



events corresponds to the cross section and acceptance 
at M t = 165 GeV/c 2 . Of these Monte Carlo experi- 
ments, 31% yielded a statistical uncertainty less than 3.9 
GeV/c 2 . 

A summary of systematic uncertainties in this mea- 
surement is shown in Table HT] The largest source of 
systematic uncertainty in our measurement is due to un- 
certainty in the jet energy scale 28], which we estimate 
at 3.5 GeV/c 2 by varying the scale within its uncer- 
tainty, including effects of high instantaneous luminos- 
ity (which have been found to contribute an uncertainty 
of 0.2 GeV/c 2 ). This is necessarily larger than in the 
previous application of this method , as we have in- 
cluded additional jets measurements in our calculation; 
future measurements would benefit from a direct cali- 
bration of the 6-jet energy scale from Z — > bb decays. 
We estimate the uncertainty due to the limited number 



of background events available for Monte Carlo exper- 
iments to be 0.7 GeV/c 2 . Uncertainties due to PDFs 
are estimated using different PDF sets (CTEQ5L [22| vs. 
MRST72 [2{|), different values of A QCD and varying the 
eigenvectors of the CTEQ6M [22| set; the quadrature 
sum of these uncertainties is 0.8 GeV/c 2 . Uncertainty 
due to showering model in the Monte Carlo generator 
used for tt events is estimated as the difference in the ex- 
tracted top quark mass from pythia events and herwig 
events and amounts to 0.9 GeV/c 2 . We estimate the un- 
certainty coming from modeling of the two largest sources 
of background, Z/"f* and events with a misidentified lep- 
ton, to be 0.2 GeV/c 2 . Uncertainty due to imperfect 
modeling of initial-state (ISR) and final state (FSR) QCD 
radiation is estimated by varying the amounts of ISR and 
FSR in simulated events 0, giving 0.3 GeV/c 2 for FSR 
and 0.3 GeV/c 2 for ISR. The uncertainty in the mass 
due to uncertainties in the response correction shown in 
Fig. [T] is 0.6 GeV/c 2 . The contribution from uncertain- 
ties in background composition is estimated by varying 
the background estimates from Table [J within their un- 
certainties and amounts to 0.7 GeV/c 2 . The uncertainty 
in the lepton energy scale contributes an uncertainty of 
0.1 GeV/c 2 to our measurement. Adding all of these 
contributions together in quadrature yields a total sys- 
tematic uncertainty of 3.9 GeV/c 2 . 

In summary, we have presented a new measurement 
of the top quark mass in the dilcpton channel, M t = 
164.5 ± 3.9(stat.) ± 3.9(syst.) GeV/c 2 . This is the most 
precise measurement of M t in this channel with an ap- 
proximately 35% improvement in statistical precision 
over the previous best measurement [3] • The systematic 
uncertainly, while 15% larger, is nearly completely corre- 
lated with systematic uncertainties in measurements in 
other channels and so does not impact the global combi- 
nation nor an analysis of measurements in different chan- 
nels. Previous measurements yielded smaller values of 
Alt in the dilepton channel 0, Q than in the single 
lepton 0] and all-hadronic [3l[ decay channels, though 
the discrepancy was not statistically significant. Our 
measurement continues that trend with substantially in- 
creased statistical precision. A global combination [32j |. 
however, shows that these variations are consistent with 
statistical fluctuations. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundesministerium 
fiir Bildung und Forschung, Germany; the Korean Sci- 
ence and Engineering Foundation and the Korean Re- 



7 



search Foundation; the Particle Physics and Astronomy 
Research Council and the Royal Society UK; the Institut 
National de Physique Nucleaire et Physique des Partic- 
ules/CNRS; the Russian Foundation for Basic Research; 
the Comision Intcrministerial de Ciencia y Tecnologfa, 
Spain; the European Community's Human Potential Pro- 
gramme under contract HPRN-CT-2002-00292; and the 
Academy of Finland. 



[1] The LEP Collaboration, The LEP Electroweak Working 
group, The SLD Electroweak and Heavy Flavor Groups, 
Tech. Rep. CERN-PH-EP/2005-051, CERN (2005). 

[2] S. Heinemeyer et al, J. High Energy Phys. 0309, 075 

(2003) . 

[3] G. L. Kane and S. Mrenna, Phys. Rev. Lett. 77, 3502 
(1996). 

[4] A. Abulencia et al. (CDF Collaboration), Phys. Rev. 

Lett. 96, 152002 (2006). 
[5] F. Abe et al (CDF Collaboration), Phys. Rev. Lett. 82, 

271 (1999). 

[6] B. Abbot et al. (D0 Collaboration), Phys. Rev. Lett. 80, 
2063 (1998). 

[7] A. Abulencia et al. (CDF Collaboration), Physical Re- 
view Letters 96, 022004 (2006). 

[8] T. Affolder et al. (CDF Collaboration), Phys. Rev. D 63, 
032003 (2001). 

[9] V. Abazov et al. (D0 Collaboration), Nature 429, 638 

(2004) . 

[10] V. Abazov et al. (D0 Collaboration), PRD (2006). 

[11] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 

73, 112006 (2006). 

[12] K. Kondo, J. Phys. Soc. Jpn. 57, 4126 (1988). 

[13] G. Goldstein et al, Phys. Rev. D 47, 967 (1993). 

[14] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 

74, 032009 (2006). 



[15] B. Jayatilaka, Ph.D. thesis, University of Michigan 
(2006). 

[16] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
032001 (2005). 

[17] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 
93, 142001 (2004). 

[18] Missing transverse energy, Et, is defined as the magni- 
tude of the vector, — J]V E T rti, where E T are the magni- 
tudes of transverse energy contained in each calorimeter 
tower i, and Hi is the unit vector from the interaction 
vertex to the tower in the transverse (x, y) plane. 

[19] The unclustered transverse energy in an event is the total 
transverse energy in the event that is measured in the 
calorimeter but not clustered into a lepton or jet. 

[20] G. Corcella et al, J. High Energy Phys. 01, 010 (2001). 

[21] T. Affolder et al, Nucl. Instrum. Methods A 447, 1 

(2000) . 

[22] J. Pumplin et al, J. High Energy Phys. 0207, 012 (2002). 

[23] G. Mahlon and S. Parke, Phys. Lett. B 411, 173 (1997). 

[24] G. Mahlon and S. Parke, Phys. Rev. D 55, 7249 (1997). 

[25] While up to 15% of it pairs at the Tevatron are pro- 
duced by gluon-gluon fusion (gg — » it), this term can be 
excluded from the matrix element with negligible effect 
on the precision of the measurement. 

[26] M. L. Mangano et al, J. High Energy Phys. 0307, 001 
(2003). 

[27] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 

(2001) . 

[28] A. Bhatti et al, Nucl. Instrum. Methods A 566, 375 
(2006). 

[29] A. D. Martin et al, Phys. Lett. B 356, 89 (1995). 

[30] A. Abulencia et al. (CDF Collaboration), Phys. Rev. D 

73, 32003 (2006). 
[31] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 

(Submitted), hep-ex/0612026. 
[32] The CDF Collaboration, The D0 Collaboration and 

Tevatron Electroweak Working Group, Tech. Rep. 

FERMILAB-TM-2355-E, Fermilab (2006). 



